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^vj We report on recent work concerning the effect which the change in vacuum struc- 

_^ ^ ture (negative energy Dirac sea), in the presence of a confining scalar field, has on 

— the nucleon structure functions and parton distributions. Using the Dirac equation 

in 1+1 dimensions, we show that distortions in the Dirac sea are responsible for 
part of the violation of the Gottfried sum rule — i.e., part of the fiavor asymmetry 
in the proton sea. Our basic argument is that, even if isospin is an exact symmetry, 
the presence of a confining potential changes the vacuum structure, and inevitably 
leads to a violation of SU (2) fiavour symmetry in a hadron with a different number 
fH , of valence u and d quarks. The same mechanism also leads to a prediction for Am 



and Ad. 



Oh! 

1) ■ 1 Introduction 



The study of the violation of the Gottfried sum ruleO has attracted considerable 
interestpSince its experimental confirmation by the New Muon Collaboration 
(NMC)H. At present, there are several possible explanations for it, including 



2auli exclusion principle at 



mesonic (mainly pionicEl) contributions and .-f|,-_,_|«^ 
the quark levellj - for recent reviews^ee Refs.EiD'LlH. 

Here we report our recent workH, on this problem, aimed at putting the 
phenomenon into a more general context. Our approach is based on a totally 
nonperturbative effect, namely the change in the vacuum structure (including 
distortion of the negative energy Dirac sea) in the presence of a confining 
potential. Such changes in vacuum structure have been studied iiL|gijea|t detail 
in other contexts (including the phenomenon of fractional chargell2lO't3), since 
the initial work of Jackiw and RebbilliJ. 



"Presented at the Joint Workshop of the Special Research Center for the Subatomic Structure 
of Matter (CSSM) and the National Institute for Theoretical Physics on Lepton Scattering, 
Hadrons and QCD, March 26 - April 6, 2001, Adelaide, Australia. 

1 



In the case of the nucleon structure functions, self-consistent studifiS-nLthfi 
vacuum structure have been made within a chiral quark soliton modelll3'ot3ii3. 
However, it is not easy to see how the change in vacuum structure is directly 
reflected in the calculated results. 

In this report we will explicitly show that the existence of a (quark con- 
fining) potential distorts the negative energy Dirac sea and inevitably leads 
to the violation of the Gottfried sum rule. We emphasize that our argument 
holds even if charge symmetry and flavor SU(2) symmetry are assumed to be 
valid. We also predict Au in proton based on the same mechanism. 

2 The Dirac equation in 1+1 dimensions 

Consider massless (m = 0) fermions in 1 + 1 dimensions, in the presence of 
a static scalar potential V{z). The fcrmion spectrum is determined from the 
Dirac Hamiltonian, 



Hipiz) = 



-ia— + pV{z) 



viz), 

ifiiz) ^ Eifiiz), (1) 



dz 

-^,+V{z) 

.ii + n^) 

where we have chosen the Dirac matrix representation with a — a2 and (3 = ai. 
We choose V{z) to be an odd function, with a kink-like profile such that 
y(oo) = Vo = — y(— oo). Then the spectrum is necessarily symmetric about 
E = 0, and there are continuum thresholds at ii^ = zLVq. There is a bound 
state at E' == 0, but it is irrelevant for our discussion. Instead, we consider 
bound states in the gap |-B|(7^ 0) < |Vo|, which occur in pairs symmetrically 
placed about E = 0. 

We choose the profile for the scalar potential: V{z) = Vo tanhz, where 
Vb > 0. Then the E ^ bound state wave functions with energy £"„ = 
±y^n{2Vo — n), are given by 



■P^°-"(tanhz) 



X^iz) = Nj^^/l '-->\ (2) 

for n — 1,2, .. . [Vb], where [Vb] means the greatest integer less than Vq. In 
Eq. (0), Pj^{z) is the associated Legendre function. If Vb is an integer then there 
are real (but non-normalizable) threshold states obtained by setting n = Vb in 
Eq. (y). The continuum state wave functions are given by 



V\ 



V ±P^_i(tanhz) ; 



with energy Ep ' — ± \/Vq + p^ . For simplicity, we choose Vb = | , where there 
is a zero energy state, and two other bound states with energies £{ — ±a/2: 

(±1/ X /\/2tanh2:sech^"2;\ /.-, 

Ai(^)-[ ±,echV^, )■ (4) 

These solutions are to be compared to the free Dirac equation solutions of 
positive and negative energies {Ep = ±|p|): 

uA^) = ^ (-^f ^") e-, .,(z) = -L (^^_/H) e--. (5) 

Then, the field operator can be expanded in terms of either {up, Vp, — cx3 < 

p < oo} or {xbmr^p I'-Pp , — oo < p < oo}, each of which is a complete, 
orthonormal set 113: 

/•oo 7 

i'(^)=j 7r:[bpUp{z) + dlvp{z)\ (6) 



27r 



dp 

enXb,n + I -Z- 



ap^(+)(z) + ct^(-)(z)l. (7) 



The creation and annihilation operators for the set with Vq 7^ satisfy, 

{eL em} = Sn„i, {al, Uq} = 2tt6{p ~ q), {c^„ Cq} = 2tt6{p - q), (8) 

and all other anticomutators vanish. The operator sets, e.g., {e„,Op,cj,} and 
{bp, djj} can be related using the orthogonality (Bogohubov transformation): 

Sn = {Xb,n\Up)bp + {Xb,n\Vp)dl, (9) 

ak = {ip^^'^\up)bp + {ipi+^\vp)dl (10) 

cl^{ipi7^\up)bp + {ipi-^\vp)dl (11) 

where the bra-ket notation implies, /^ ■^. 

Let us focus on the bound state, \xb,i), where the wave function is nor- 
malized to unity and its expression is given by Eq. (0). Using the relation of 
Eq. (H) and the property of orthonormality, the normalization of the bound 
state can be written by: 

1 - (xmIxm) = r ? (I(xmI"p>I' + \{Xb,i\vp)n . (12) 



In Eq. (p^), the term |(x6,ikp)P i\{Xb,i\up)\'^), can be interpreted as the 
probabihty density distributions with which the bound state |xb,i) would be 
filled by the negative (positive) energy free states, if the potential were turned 
off suddenly (but not adiabatically). 

We show in Fig. |l] the quantity |(xft,i|"p)P + l(x&,il^p)P' (square of the 
wave function overlaps) as a function of energy for the free plane wave solu- 
tions, Ep(free). Surprisingly, one can immediately notice that the positive 
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Figure 1 : Bound state number density distributions as a function of free plane wave energy, 
Ep{free), when projected onto free plane wave solutions. 



energy bound state wave function, Xb,i, has non-zero overlaps with the nega- 
tive energy free plane wave solutions (the shaded area in Fig. |l|). If we integrate 
over Ep{free), it gives unity by construction, or probability conservation, as 
it should be. We should also mention that the negative (positive) energy con- 
tinuum state wave functions of Eq. (j^ have also non-zero overlaps with the 
positive (negative) energy free plane wave solutions. 

The conclusion of this section is that the positive energy bound state, 
appearing in the Fock space built on the presence of the potential, has non- 
zero overlaps with the negative energy free plane wave solutions. It is this 
point that we want to make a connection with the deep-inelastic scattering 
(DIS) and the parton model. 



3 Connection with the Parton Model 

In this section we consider the interpretation of the results obtained in the 
previous section in connection with the parton model. It may be useful to 
recall th£ assumptions which leads to the standard parton model according 
to Ref. tZI: (1) The current couples to quarks. Then the contributions to the 
forward Compton amplitude can be classified by the flow and interactions of 
quark lines. (2) At large value of Q^ the currents, hut not the states, may 
he treated as in free field theory. Thus, final state interactions and vertex 
corrections are ignored. 

Let us now identify the positive energy bound state wave function of Eq. (^ 
as that of the valence quark in a simple, model proton. We recall that in 
Fig. |l| the quark number density distribution is projected onto free plane wave 
solutions of energy Ep(free). As mentioned above, these free plane wave 
solutions are now to be identified with the partons, which appear in free field 
theory. As expected, the valence quark wave functions in the proton have 
non-zero overlap with the negative energy free plane waves, and this means 
that the distributions are non-zero in the region where Bjorken x is negative. 
As a result the integral of the valence distribution over the physical region 
X > will not be one. The resolution of this problem lies in the fact that 
the Dirac sea in the presence of a scalar potential is also different from that 
of free space. In particular, the Dirac sea in the presence of the soliton in our 
ID example (or the confining potential for the proton) consists of having all 
negative energy states occupied. If re-expressed in terms of free plane waves, 
however, one finds an equal number of positive energy free states occupied 
and negative energy states (holes) empty. That is, when the occupied negative 
energy states are expressed in terms of free states one finds an intrinsic sea of 
qq pairs. In the absence of valence quarks this sea carries no net flavour (or in 
3D, spin), relative to the free space vacuum. 

However, once we add a single valence quark (say a u quark), the negative 
energy tail will fill some of the holes in the (free space) Dirac sea. An equal 
number of occupied positive energy states from the u quark Dirac sea will then 
be interpreted as the missing valence normalization (because Uy is defined to 
be u — u). The net result is that we have one valence u quark, together with 
a reduction in the number of uu pairs in the sea of the proton. That is, one 
will have d > u by an amount equal to the area under the valence curve to the 
left oi Ep = in Fig. 1. If instead we had one valence u and one valence d the 
intrinsic sea would be flavour symmetric again, while for two valence u's and 
one valence d, as in the proton, we have d > u hy the same amount as for the 
single valence u (the area under the valence curve to the left oi Ep = in Fig. 



1.4 





C\J 


1/. Dnlp-nrp. 


^ 1.2 


■ d- valence \ '■ 


i 


\ 


Q, 1.0 


■ \ - 


CD 


\ 


s>. 


\ 


^ 0.8 


1 \ 


^ 




\ 


•+-3 






1 0.6 


, 




fi, 


1 


/ \ \ 


ti 0.4 


- 


/ \ \ 


5- 




\ \ 


Q) 


f 


\ v 


g 0.2 


1 


\ \. 


o 


1 


\\ 




# 


\ ^^^^ 


0.0 


: _-m«in— if??*!^ /. 





Ep(free) (dimensionless) 

Figure 2: u and d valence quark number density distributions in the "proton" as a function 
of free energy, Ep(free), projected onto free plane waves. 



!)■ 

The situation with respect to the Gottfried sum rule is therefore iUustrated 
in Fig. y. Because there are more valence u's, the u sea is depleted more than 
the sea of d's. Mathematically, the expression for the Gottfried sum rule is 
then: 



[Frw-Frw], 



1 dx 

^ 

1 /•! 2 [^ - 

- dx[uy{x) - dy{x)] - - dx[d{x) - u{x)], 
-J Jo •J Jo 

1 2 

3^3 



dp 

2^ 



\{Vp\XbA)? 



(13) 



(Note that we have assumed charge symmetry to set d" = u^ = u, etc.) 

4 Conclusion 

By considering the simple ID problem outlined here, we aimed to make the 
physics origins of one source of flavour asymmetry in the proton explicit and 
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easy to understand. While one cannot access spin in the ID case, the simpUcity 
of the argument for flavour makes the corresponding argument for Am (and 
Ad) clear. Since in a spin up proton described by an SU(6) wave function 
the valence u quark is predominantly spin up, one will be missing spin up u 
quarks in the sea of the proton. Thus Am will be jpnsitive and of a similar 
magnitude to that of d — u from the same source cl'Ea. Since the d quarks 
tend to have spin down in the proton, one also expects Ad to be negative and 
considerably smaller in magnitude than Am. Since any contribution to 
observed flavour asymmetry from the nucleon's pion cloud is unpolarizedt 
the experimental determination of Am and Ad would be extremely valuable in 
reaching a quantitative understanding of the roles of these two very different 
physical mechanisms in generating the non-perturbative structure of the proton 
sea. 
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